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ABSTRACT 


Hypervelocity flows for velocities in excess of 1.4 km/sec (Mach 5) require 
very high stagnation temperature to avoid liquefaction. The arc heater wind tunnel 
has been designed to provide such flows. The electric-arc driven wind tunnel can 
develop stagnation temperatures up to 13,000°K which will produce hypervelocity 
flows up to 7 km/sec (earth orbital speed). The nature of the flow, however, is such 
that the high temperature source flow may cause severe gradients at the nozzle exit. 
In order to perform aero-thermodynamic tests the characterization of the flow in the 
test section is required. This paper experimentally determines the stream profiles for 
an arc jet wind tunnel conical nozzle directly from calorimetry and pitot probe 


Surveys. 
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I. INTRODUCTION 


This report is concerned with developing a procedure by which flow field 
characteristics of very high temperature arcjet flows may be produced directly from 
stream surveys of nozzles. 

Hypervelocity flight in the atmosphere has previously been concerned with 
entry into the atmosphere and descent to the ground. This has been expanded 
during the last decade to include vehicles which will use the upper layers of the 
atmosphere to perform hypervelocity aerodynamic maneuvers and vehicles which 
will fly from the ground into orbit using air breathing powered lift. The design and 
testing of an aeroassisted orbital transfer vehicle (AOTV) and the National Aerospace 
Plane (NASP) will require ground based testing in order to obtain a data base for 
this flight regime. Such data will be used to validate 3-D, reacting real-gas, CFD 
codes that will be used in the design of advanced vehicles where aero- 
thermodynamic effects are important. 

The 60-MW Shuttle Interaction Heating Facility (IHF) was designed to meet the 
testing requirements for development of the thermal protection system of the shuttle 
orbiter vehicle. Its design provided for heat transfer simulation and stagnation- 
region flow simulation for conditions which the shuttle orbiter would encounter. If 
this facility is to be used in aerodynamic testing, the arcjet flow characteristics need to 
be known in order to determine if they provide optimum simulation of the actual 
flow conditions. 

In arcjet operations, extremely high temperatures are generated. For example, 
for the test reported here, the stagnation temperatures are of the order of 6000°K 


(10,800°R). Probes must be highly cooled in order to survive the environment. 


Conventional techniques using thermocouples, for example, to determine local 
stream energy simply will not work because of these high temperatures. Stream 
tube flow-swallowing techniques can be used up to moderate temperatures, 
however, such enthalpy probes require sharp leading edges; and these can not be 
cooled in the arcjet environment. Moreover, such probes with blunt cooled leading 
edges yield an uncertain stream-tube capture; and resolution and accuracy both have 
large uncertainties. To circumvent these problems for arcjet stream surveys, a simple 
method has been developed that uses small diameter, cooled, hemispherical-shaped 
probes to measure local heat transfer and pitot pressure. It is shown in this report 
that these two measurements yield stream profiles to a high degree of resolution. 
The theoretical method derived here has been formulated into a computer code for 
data reduction of stream surveys. This code can be incorporated into a real-time data 
acquisition system and will provide on-line reduced stream profiles during the 
nozzle survey for immediate assessment. 

PROFILE is a computer code developed to use the 60-MW IHF stream survey 
data for a 13-inch conical nozzle. The PROFILE code directly produces normalized 
velocity, density, mass flux and enthalpy profiles. The derivation of PROFILE uses 
the basic conservation equations and is described in Chapter IV. PROFILE uses the 
heat transfer and stagnation pressure data as inputs and outputs normalized values of 
heat transfer rates, stagnation pressure, velocity, density, mass flux and enthalpy 
with respect to centerline values of the flow. Characteristics of flows for other arcjet 
designs and other nozzles can be readily obtained using heat transfer and stagnation 


pressure survey data in conjunction with the PROFILE code developed here. 


II. INTERACTION HEATING FACILITY DESCRIPTION 


A schematic diagram and photograph of the IHF arcjet facility are shown in 
Figure 1 and 2, respectively. The major components of the facility include a 
constricted-arc air heater, a 150 MW power supply, two interchangeable hypersonic 
Mach number nozzles, a walk-in test chamber, and ancillary subsystems consisting 
of a steam-ejector vacuum system, cooling water system, and automated data 
acquisition system. The high-power capability of the arc heater (rated to 60 MW) 
provides for full scale subsystem tests in high enthalpy streams in both large scale 
stagnation flows and boundary layer flows at high Reynolds numbers. 

The primary component of the facility is a constricted-arc heater 390 cm 
(153.3 in.) long and 8 cm (3.15 in.) in diameter shown in Figure 3. The rated 
airflow rate is 1.4 kg/s (3.09 lbs/s) and the rated pressure is 12 atm. Cooling water 
is supplied at 7000 kPa (1000 psi), which provides cooling for a maximum heating 
of 12 kW/cm® to the constrictor wall. Eight anode and cathode electrodes, each 
capable of conducting 1000 A, provide a maximum arc current capacity of 8000 A. 
Thus, the capability for continuous operation at high current and high enthalpy 
levels is provided. 

Direct current for operating the arc heater 1s provided by a 150 MW D.C. 
power supply. Six modules make up the power supply, each consisting of a three- 
phase, full-wave, phase controlled silicon rectifier. One module is capable of 
operating up to a level of 5500 V into an arc load at current levels up to 2700 A. 
The six modules can be connected in any compatible series/parallel arrangement by 


means of remotely operated setup switches. The duty cycle for the power supply is 


30 minutes on, 30 minutes off. A 150 MW rating is achieved for run durations of 
approximately 15 seconds. 

The facility has two interchangeable nozzles: 1) a semielliptic, Mach number 5 
nozzle with an open-jet, flat-plate test section measuring 80 X 80 cm (31.5 X 31.5 
in.) for boundary-layer surface flow tests and 2) a conical, nozzle for testing large- 
scale stagnation flows in a free jet. The conical nozzle is constructed in 5 sections 
with exit diameters varying from 15.24 cm (6 in.) to 105 cm (41.4 in.). The Mach 
number varies with the nozzle area ratio from 3 to 7.5 depending on the exit 
diameter selected. Both nozzles are water cooled by internal water passages, and the 
throat diameters are 6.03 cm (2.37 in.). Both nozzles are designed to flow into an 
evacuated, walk-in test chamber that can be maintained at a pressure of less than 
0.004 atm (3 mmHg), a pressure that varies with flow rate according to the 
pumping curve of the steam ejector system. Tests have shown that supersonic flow 
with 33 % blockage ratio has been attained in the conical nozzle, corresponding to a 
model diameter of about 60 cm (24 in.). 

An on-line automated data acquisition and computing system is used to 
record arcjet performance and test data. The conversion from analog to digital 
signals is made by means of an integrating digital voltmeter, so that the signals can be 
processed by an on-line computer. Output consists of a real-time display of five 
channels that are monitored during a run, printed and reduced data, and a 


continuous time history data plot. [Ref. 1] 


II. SURVEY DESCRIPTION 


The design of the IHF allows for the control of two variables: 1) constrictor 
pressure and 2) arc current. The test conditions are uniquely determined by these 
two control settings only. For this survey, the test chamber pressure was 
maintained at a low pressure to allow the measurements to be made in a free jet. No 
attempt was made to control the expansion ratio into the free-jet chamber. In general, 
the flow was as an under-expanded free jet. This type of operation causes the free jet 
to spread and to operate at a slightly higher Mach number relative to the nozzle exit 
plane.The survey measurements close to the nozzle exit plane, then are relatively 
uninfluenced by the free-jet expansion process. Most of the discussion pertains to 
the 7.62-cm (3-inch) station to avoid the complexity of the additional expansion. 
Constrictor pressure was varied for several values of current in order to map the 
operating envelope for stagnation flows using the conical nozzle. 

The probes used to obtain data consisted of a hemispherical Gardon type 
calorimeter and a pitot probe shown in Figures 5 and 6. The probes were 1.58 cm 
(0.625 in) in diameter and were mounted 20.32 cm (8.0 in) apart horizontally. The 
probes were mounted to a transverse mechanism that was positioned to pass 
horizontally through the flow. Sweep measurements were taken at 7.62 cm (3.0 
in), 19.5 cm (7.7 in), 34.8 cm (13.7 in) and 52.6 cm (20.7 in) downstream of the 
nozzle exit. 

Figure 7 shows probes, probe supports and their relationship to the nozzle 
during a test using a 15.24 cm (6.0 in.) exit diameter nozzle. Figure 8 shows 


probes entering the flow. Figure 9 shows probes approximately on centerline of 


nozzle. All photographs were taken using only radiated light from gas behindthe 
normal shocks on probes and support mechanisms. 

Data consisted of heat transfer rate and stagnation pressure measurements taken 
at 0.212 cm (0.083 in.) intervals across the flow. Typical profiles for heat transfer 


rate and stagnation pressure for each current setting are shown in Figure 10. 


IV. DERIVATION OF DIRECT METHOD FOR PROFILES 


Determining the stream profiles for an arcjet wind tunnel nozzle requires that the 
heat transfer rate and the stagnation pressure measurements behind a normal shock 
make up the only parameters on which the velocity characteristics depend. Briefly, 
this method is based on the fact that stagnation-point heat transfer is proportional to 
the product pV3 (approximately) while the pitot pressure is proportional to pV2. 
When these two quantities are known, a simultaneous solution gives density and 
velocity directly from the measurements. Once the velocity profile has been 
determined the profiles for enthalpy, mass flux and density can be readily found. 

In the following derivation, free stream flow parameters are those found in front 
of the normal shock and will be designated by subscript one while subscript two 
will denote parameters behind the normal shock. The starting point is the 
momentum and continuity equations in the following form: 

P,+ PV, = P,+ PW, (1) 
pv, =PY, (2) 


Rearranging the momentum equation and substituting into it the continuity equation 
] _ la ee ee ee 
nf ole 
Fa Py (3) 


2(1-(1/(P,/P ,)) 
Pa Pa ; 7 (V/( Pp,’ i 


gives: 


Solving for p2 


Noting that M2 is much less than 1, that M2 = @( y, M) and recalling that from state 2 


to Stagnation state 2 is along an isentrope then: 


Ke a (1+ (4 *)mz) (5) 





Substituting equation (6) into a normalized equation (4) gives: 


y? 
ie PM 1 


("%) (pW on 


. (7) 
At this point the normalized stagnation pressure behind the shock is related to 
the velocity in front of the shock. The next step introduces the closed form of the 
Fay-Riddell [Ref 2] equation derived in Appendix A., which is a correlation for the 


heat transfer rate to a hemisphere; namely, 


P; 
'= 00417. /—*(h, - 
: R (%~ %) (Btu/sec-ft?) (8) 
where hy is the total enthalpy of the flow and hy is the enthalpy of the wall. The 


units are: h, Btu/lbm; py, atm; and R, ft. Normalizing with centerline values gives: 








(9) 
CL 
(9) 
Using the following form of the energy equation 
Ve Vv; 
a er ee 
hts rhtz=h, 10) 


and equation (7) in equation (9) while neglecting the enthalpy of the wall compared 
to the total enthalpy of the flow yields: 


pV (V" + 2h) 























ee 
(q) Vv? Vi 4+2h 
CL (p, om ( 1 a > CL (11) 
Rearranging 
2h, 
3 La 
q _ Pp, ve | ( Vi 
ce Poe, = (1+ 7”) 
Ma 
CL (12) 


Since the velocity in front of the shock is supersonic the last term drops out. 


Changing notation to reflect only centerline and free stream gives: 








(13) 


Squaring equation (13) and dividing by equation (7) removes the density term and 


results in an equation which relates free stream velocity, stagnation pressure and heat 


4 omnes =) 
(7) = (ez) (FS) 
Vor Ici) \ Pt CL (14) 


This equation equates the normalized velocity of the free stream to the two measured 


flux. 





quantities behind the normal shock. 

The assumption made to derive this result are simply that the flow is 
axisymmetric and the enthalpy of the wall is small compared with the total enthalpy 
of the free flow. The wall enthalpy is measured directly from cooling water and is 
found to be approximately 0.3 MJ/kg (134 Btu/lbm) or less than 5 % of total 
enthalpy. 

Once the normalized velocity has been found it 1s used in equation (7) to find the 


value for the normalized density. 





(15) 
Normalized mass flux is obtained from the product of normalized density and 
velocity. 


sii amg (F—} 


Mier Per ee 


(16) 
The normalized enthalpy is taken directly from the energy equation where 
gravitational effects are neglected and the static enthalpy is assumed to be much less 


than total enthalpy, which results in the following equation. 
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hi (~_) 
Ce Von (17) 
Eguations (14), (15), (16), and (17) were used in the code PROFILE to calculate the 


characteristic profiles of the arcjet high enthalpy flows using survey data. 
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V. RESULTS 


For high-powered, arc-heated, hypersonic facilities the testing for aerodynamic 
or thermodynamic (heat transfer) studies requires a central core flow with properties 
well characterized. Ideally, of course, constant core flow properties are desirable. 
The present study was undertaken to characterize the flow properties of the 60-MW 
arc heated wind tunnel at NASA-Ames Research Center. Such flow properties are 
essential to the engineer for interpretation of the results for aero-thermodynamic 
tests. 

Arc heaters operate with a large enthalpy gradients near the walls. These 
gradients cause the flow properties in the central core to vary; and are an important 
concern in simulation in these ground-based test facilities. The degree of this 
variation has not been reported up to the present time for the arcjet facility at NASA- 
Ames Research Center. Enthalpy gradients are inherent in all types of arc heaters 
that operate with cooled walls. For example, the temperature ratio of the wall for a 
typical arc heater varies from 2 - 5 % when the enthalpy level is in the typical vehicle 
reentry range of 18.6 - 41.87 MJ/kg (8000-18000 Btu/lbm) which in terms of 
velocity equates to 20,000-30,000 ft/sec. The constricted arc heater was developed 
to minimize this gradient by the current density and length to cause the enthalpy 
profile to assume a flatter shape than a free-burning or vortex-stabilized arc heater. 
The theory of Stine-Watson [Ref. 3] shows that for a constricted arc heater the 
enthalpy profile assumes a Bessel function form as opposed to a linear profile for a 
free-burning or vortex-stabilized design (Huels arc) illustrated in Figure 11. 

Results from the profile measurements obtained here can be used to answer the 


question as to whether the property variables (enthalpy, velocity, density and 
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energy flux) have the same profile or shape at the nozzle exit as they do for the 
constrictor. Scant measurements exist for the constrictor. Most of the values for the 
constrictor arc heaters are predicted by the Stine-Watson theory. Elaboration of this 
theory is represented by a computer code called ARCFLO I [Ref. 4]. This code 
was developed to predict the performance of high-pressure, high-enthalpy 
constricted-arc heaters. The conservation equations (mass, axial momentum, and 
energy) are written for a steady-state and radiatively participating gas mixture with 
electric current conduction. The governing equations employ the boundary-layer 
approximation and, hence, the pressure in the wall normal direction is assumed 
constant. Auxiliary relations, such as the equation of state, global current continuity, 
and global mass continuity, are used to relate the gas density with the pressure, mass 
flow rate, and mixture properties. 

The profiles of density, velocity, enthalpy and energy flux computed by 
ARCFLO I in the constrictor are compared with experimental survey data obtained 
at the nozzle exit in Figure 12. The difference seen in the density profiles, Figure 
12a, can be accounted for by the increase in density near the wall of the constrictor as 
opposed to the decrease in density in an expanding free jet. The velocity profiles, 
Figure 12b, are very similar with the experimental data showing a slightly flatter 
profile. The enthalpy profiles, Figure 12c, shows for ARCFLO II a central core 
region where the enthalpy drops moderately from the centerline to r/R = 0.4, a 
linearly decreasing region to r/R = 0.95, and then a steep drop near the wall where 
the conduction losses dominate. The ratio of centerline enthalpy to mass average 
enthalpy is 1.6. This low ratio value is typical for wall-stabilized arcs. The survey 
data exhibits a much flatter profile than ARCFLO II with a ratio of centerline 


enthalpy to mass average enthalpy of 1.15. 
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The product of the above three profiles is shown in Figure 12d and is the 
energy flux of the airstream. Note that this parameter is proportional to the 
convective heat-transfer rate to a body within the stream and that it is a measure of the 
usable stream area. The flat profile for ARCFLO II implies that stagnation point 
heating to a test body would be constant to within 10 % of the value obtained on the 
centerline flow for r/R up to 0.7. The survey data profile shows a more 
pronounced curve which can be attributed to the decrease in density in the free jet. If 
test chamber Static pressure is maintained at the nozzle exit static pressure then the 
density profile would become flatter and thereby improve the energy flux profile. 

Figures 13 - 28 are plots of PROFILE code output data at specified survey test 
points. All test point plots are normalized with respect to the corresponding values 


obtained at the 3-inch downstream positions. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The present study of survey data for the 13-inch IHF nozzle showed: 

: That there is good agreement between the ARCFLO II code 
prediction for velocity profile and experimental data; 

- Enthalpy profile of experimental data is much flatter than ARCFLO I 
theory predicts, which can be explained by the codes inability to 
account for energy transfer toward constrictor edges. 

It is recommended that: 

- Future surveys control test chamber conditions in order to determine 
if improvements in the profile characteristics result; 

- Sting mounted heat transfer and stagnation pressure probes be 
developed in order to reduce probe mount interference with 
boundary layer of the flow; 

- Survey data be organized into accessible files so that characteristic 
maps for individual nozzles and complete nozzles can be utilized by 


researchers and contractors. 
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APPENDIX A 


DERIVATION OF THE CLOSED FORM OF THE FAY-RIDDELL EQUATION 
by WARREN WINOVICH 


1. The Nusselt number is defined by: 





Nu = (w = wall value) 
K, (Al) 








H(h, — h,) 
¢ = HiGe— te = ae 
r (A2) 
where 
— Cp rec Se h,-h, =] 
ais Cp “Pp 7 Lalas aa 


(by definition of Nusselt number; wall temperature used as reference) 


3. Solving for heating rate in terms of Nusselt number. 


Nuk, (h, -—h,) 
7 = nace 


Pe (A3a) 


, 1 Chas 
at ese ~ VP wb (2 C, xy 
Hw k, (A3b) 


where subscript ‘e’ indicates external flow conditions 
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j= eI nS 











(A3c) 

4. For the laminar boundary layer near the stagnation point: 

€ fed 
u, = x 
dx, (A4a) 
therefore 
Sa pe iP a 

7 Pw Hw 5 Cc 

(A4b) 


5. Evaluation of velocity gradient at the stagnation point for a hemisphere: 





=p (momentum equation ) 








BePuy  P. 
u, = 5 t — > 
e JA (ASa) 





P. : 
= 1-sin?(=) (Newtonian flow over hemisphere ) 


= 4 ( =] sin(z) (ASb) 
a mn 22, 5+) eos(2) (ASc) 



































P : (A5d) 
6. The heating rate equation becomes: 
(5+ | (h, = (h, — hb.) 
Pre, (A6) 
7. The Fay and Riddell result for the heating transfer coefficient 1s: 
N 0.4( Dalle \ 
u ‘ sks 
= 0 .763(Pr) ( c 
V/ Re Pp. H. (A7) 
(equation 63 and 45 Fay and Riddell) 
8. Substituting equation (A7) into equation(A6); using 
bare 
- wR Dean and R=535 
4/\ P, 
gives: 
G = 0.763( Fri” 4 2eS Gna ONL. 
0.4 
oon “rem a0) 
poi "VR (A8a) 
qv R 0.763 
Ooty = D2 ote om (GEE) IZ a 
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9. Using: 


P, 
w= 39.65 0 
P T,Z, (A9a) 
@.628 
-7 3 
Bp = Heres (saga) cat (A9b) 
where 
LU exact 
C= an 
Hiectieils 


Equation (A9b) is appropriate for temperatures of 6000°K (10000°R) where arcjet 
operate and for reentry and reentry stagnation conditions (Ventrry > 5000 m/sec ~ 


15000 ft/sec) 


0.628 


— 39.65x184x 10 | PT, 








dni Ss 10000°°° Tez: K, 
me 0.628 
—~ -9.60x 10 ‘—-— K, 
P, L2G: (A9c) 
Py Ul. a TAZ, i’ K, 
Dal. = | a Ae 0.628 K . 
T, (A9d) 


Note that Ky = 1.00 for the expression chosen for viscosity. 


10. The final form for the expression for heating rate becomes: 


qVR — 0.763 , 2g (53.3) 3.10x eee rie” 
eh) (pry (az) 
8.1 @.0628 @.1 


DZ ees 
SD GR aan ae. 


Re) (A10a) 


hs 


Simplifying for Pr = 0.71 the equation becomes: 


6.4 


_ tva 2.) Tae 
a/p,(h,- h,) ~ (T.Z.)°" 





. @.161 6.251 

gqVR (T,Z,) Z. a4 
=~ = {constant } ———__,--> K, 
VP. (h,— hy) (T.Z,) °° \2: 


Noting that Zy = 1.00 
For Pr = 0.71 the constant in brackets above becomes: 


{constant pe 4/22 (53.3) 3.10x 10° = 0.02225 
(Pr) 


and the expression for heating rate is: 





oi 4 = (o.o2228)/x0(E2-) a 
| ee | (A10b) 


The term within the braces is only a function of the stagnation state. For usual wind 
tunnel tests the stagnation state is within the limits: 

0.01 < ps < 1 atmosphere 130 < hg < 20,000 Btu/lbm 
For this range of pressure and enthalpy the term in the braces is essentially a 


constant: 
0.0372 0.0638 


=1.873+ 0.042 
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Therefore the hemisphere formula for heating rate is: 





"5 = {0.02225}{1.873} = 0.0417 
or: 
g = 0.0417, /—=* (h,- h,) —" 
where 
g = Btu/s: ft’ 
p,=atm — p, (measured) 
R = feet | 
h= Btu/Ib,_ 
Symbols Subscripts 
Cp Specific heat e External Conditions 
h EntahIpy S Stream Conditions 
IK Thermal Conductivity t Stagnation Conditions 
Pr Prandtl Number w Wall! Conditions 
R Radius 
18 Temperature 
Zz Compressibilty Factor 
0 Density 
Ll Viscosity 
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APPENDIX B 


UNCERTAINTY ANALYSIS FOR DERIVATIONS 


Uncertainty analysis of the heat transfer and stagnation pressure probes is based 
on manufacture's claims for uncertainties. In the case of the heat transfer probe the 
manufacturer list's accuracy of readings between 3 & 4 % of the value recorded. 
For the following analysis the value of 5 % is used. The accuracy of the stagnation 
pressure probe has been found by NASA personnel to be on the order of 1/4 % of 
the recorded value. Using these values the uncertainty of the normalized velocity is 
Shown below. 


Starting with equation (14): 


F 2 “1 
Fa) “EF: 
“Gi Ger) \ Prey 


and taking the loganthmic differential gives: 











Integrating 
dv. _ 3s dg 441 -1(% pet | 
VoL 2\ q Gan nue AL Per 
dV 1 : | 1 
ae: ~ 2 [ce gam )4Icx] ~ 4 [(€r)P — (&r) Pez] 


where (er) iS uncertainty. 
In the worst case where centerline values are opposite in sign of the original 


values results in 


Je 





eV = + ((0.05) - (- 0.05)] - + ((0.0025) ~(- 0.0025)] 


dv _1 oil 
1 (0.10) - + [0.005 





so that 
GV isa 


VoL 


Applying the velocity uncertainty to the other property equations (15), (16) and 


(17) gives: 





aS 


PEPE 


AOOQE 


=H ONO 


Oo 


AAMONQANQANAAANAYNS 
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APPENDIX C 


COMPUTER PROGRAM FOR DETERMINING PROFILES 


PROF2 FORTRAN Al 


PROGRAM PROFILE 


THIS PROGRAM IS DESIGNED TO CALCULATE FLOW PROFILES FOR 

NASA AMES IHF ARCJET WIND TUNNELS. VELOCITY, DENSITY, MASS FLUX 
AND ENTHALPY PROFILES NORMALIZED TO CENTERLINE CONDITIONS 

ARE DEVELOPED. 


DIMENSION QDOTC(1000), QDOTNC1000), PTOTC1000), VNORC1000), 
& 


RHOC1000), MFLUXC1000), HTOTC1000), DISTC1000), 
QSQRC1000), PTOTNC1000) 


& 
REAL VNOR,RHO,MFLUX,HTOT,Q@SQR,PTOTN, QDOTN 


INTEGER S,1I,J,K,L,IMAX,IQDOT,IPTOT,VAL,I1,12 

REAL PMAX, QMAX,VMIN 

CHARACTER *36 G 

CHARACTER ®1 E,Y,N 

CHARACTER X*6 A,B,C 

CHARACTER *45 F 

HHH HE HEH HE HEHE HEH IE HEI HE HE HEH HEH HEH IE HEH HH HH HEH HH HE HE HE HH HHH HH HE HHH HH HH HHH HH HH HHH HK 
THIS SECTION COPIES TEST CONDITIONS INTO OUTPUT FILE 

PPS SSCS SSS SSS SS SSSSS SSCS SS SESS SS SSCS SSS SSeS SSS TS eS SSE SESS ESSE SS SS 
READ(5,1)F 

FORMATCA 


2 wy 


( 

VIRITE 
WRITEC/,1)F 

CONTINUE 
2 2. 2.2,.2.2,2,2,5.2,0,.2,0,2,%.2,0,.2,2,0,2,0,.5,5,.0,0.2,5,5.5,2.2,0,. 5,5, 5,2,5,2,5,5.5,5,5.5,0,2,5,5,5,.25,5,5,5,,5,5,5,2,5,5,5,5,5,5)5 
THIS NEXT SECTION READS TEST DATA FOR QDOT AND PTOT. 
FIRST: NORMALIZES VALUES OF QDOT AND PTOT 
SECOND: LOCATES MAXIMUM VALUE TO USE AS CENTERLINE 
THIRD: CALCULATES VALUES FOR VELOCITY, DENSITY, MASS FLUX AND 
ENTHALPY PROFILES NORMALIZED TO CENTERLINE 


NOTE: QDOT AND PTOT DATA SHOULD CONTAIN ONLY ONE ZERO AS LAST 
VALUE. THIS IS CRITERIA USED TO TERMINATE READING OF DATA. 


HEH HHH I I HE HEE I HE BE I IE IE 3 IE IE DE OE IE HEI I HE HEHE IE HE EH IE IE IE HE IE HEI TE IEE HEE HE HEHE HE HEHEHE HIE HEH HHH HE 
IMAX=0 
DO 20 I=1,1000 
READC5, 15)QDOTCI), PTOTCI), DISTCI) 
IFCQDOTCI).EQ.0)GOTO 21 
IMAX=IMAX+1 
CONTINUE 
WRITEC%,*) "NUMBER OF DATA POINTS READ =',IMAX 
ERDPIGECS 
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Oa 


AQAA 


2 
2 


2 
2 


QaDT 


G 


He HE IE IE IE EH IE IE I TEBE I HEE 3 IE IE DE I I DE IE IE IE IE I HE EI IE DE IE I IE IE IE IE DEI IE BEI I IE IE IE I HE HEHE I HE IEE I IE EK 
1 NORMALIZATION OF QDOT AND PTOT 
HEHE HE HH IH HH III HEI HEH HE HE I I IH HIE IIH I HEH HEH I IE IIH HIE HH HE I II HEH HI HH IH HEH 
QMAX = 0 
PMAX = 0 
DO 200 I = 1,IMAX 
QUAX = MAXCQMAX, QDOTC(I)) 
PHAX = MAXCPMAX,PTOTCI)) 
200 CONTINUE 
HEHE HE 3 HE IE HE IE SE 3 HE IEE I BE IE HE DE DE IE IC BE IEE DE IC HE ICE DE IE IC DE IE DE DE IE IE IE I 3 IE IE IE DE IE I IE IE IE IE IE IE IE IE I IE IE IE DE IE IE HE 
2 DETERMINING LOCATION OF MAX VALUES. THESE WILL BE USED TO 
ALIGN PROBE VALUES TO CENTERLINE. 
HE HE IE 3 IE HE HE DE 3 3 HE DE DE IE IE IE HE DE 3 3 IE IE DE SE IE DE HE SE DE 3 HE HE IE DE IE IE HE HE IE IE IE HE IE IE HEHE HE IE IE IE IE IE HE IE IE HEH HE HE HE HEHE HEE 


DO 201 I=1,1MAX 
IQDOT=I 
VAL=DIMCQMAX, QDOTCI)) 
IFCVAL.EQ.0)GOTO 202 
0) CONTINUE 


02 DO 203 J=1,1IMAX 


IPTOT=I 
VALL=DIMCPMAX, PTOTCI)) 
IFCVAL1.EQ@.0)GOTO 205 
03 CONTINUE 
05 WRITECX,*®)'IT AT QMAX =',IQDOT,' I AT PMAX =',IPTOT 
2.2 2.2.2. 2.2 2.2 005 252 £2 222 £2 ££ SF £25 SS SESS SESS SS SSS SS SSS SSS SSS SSS S SSS SS SS 
3 CALCULATION OF PROFILES 
HHH HH HK HHH HE HH HEHE HE HE HE HE HEHE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE IE HE HE HE HE HE HE HE HE HE HE HE EE HEE EEK HMR 
T1=IQDOT+179 
I2=IPTOT+179 
Pinte ce,*)*'I)] =", 1);5° I2 =',1I2 


K=0 

DO 300 J=IQDOT,11 
K=K+] 
QDOTNCKI=ABSCQDOTCI)/7QMAX) 
QSQRCK)=QDOTNCK) X¥*2 

00 : DISTC(K)=DISTCI)-DISTCIQDOT) 
=0 
DO ae I=IPTOT,I2 


=K+ 
30] PTOTNCKI=ABSCPTOTCI)/PMAX) 
DO 302 I=1,K 
VHORCI)D = SQRTCSQRTCQSQRCII/ZPTOTNCI))) 
RHOCI) = PTOTNHCID/CI/VNORCI) ) ¥X2 
MFLUXCI) = RHOCI)XVNORCI) 
S02 HTOTCI) = VNORCI)¥X2 


00 VIRITEC7,*)'QHAX = ',QMAX 
WRITEC6,*)'QHAX = ',QMAX 
HIRITEC7,%)"PMAX =',PMAX 
WIRITEC6,%)"PMAX =',PMAX 
WRITEC7,X) PssessssesesssssssssssssseltsStSssssssSsssrsrsssrssesssses 
MRITEC?7, x) * HEATING STAG VELOCITY DENSITY MFLUX 
&SENTHPY DIST! 
WRITEC7,%)' RATE PRESSURE! 
PIRITEC 7, X) 'Sastssssssssssssrsssssssssstssssssltssssssssssssssessecs 
WIRITEC7,600 
&8T(I),1=1,K 
600 FORMATC6C5 
STOP 
END 


~<-~ 
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Figure 5 Gardon type calorimeter probe. 
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Figure 6 Stagnation pressure probe - with support mount. 
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Figure 7 Probes and support mounts outside 6-inch nozzle free jet flow. 
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Figure 8 Probes entering 6-inch nozzle free jet flow during survey. 
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(a) Heat transfer rate verses radius at 3-inch 
position from nozzle exit 
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(b) Stagnation Pressure verses radius at 3-inch 
position from nozzle exit 


Figure 10 Typical profiles of the heat transfer rate 
and stagnation pressure, 13-inch nozzle (A/A* = 30) 
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Figure 11 Comparison of enthalpy profiles for 8-cm 
constricted arc. 
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— ARCFLODI 





—O— SURVEY DATA 
3 
2 
ARCFLO 0 
1 SURVEY DATA 
0 
0.0 0.2 0.4 0.6 0.8 1.0 
r/R 
(a) Density Profile (b) Velocity Profile 
SURVEY MASS AVERAGE ARCFLO LI MASS AVERAGE 
bh CENTERLINE CENTERLINE 
h 
ARCFLO D __MASSAVERAGE SURVEY MASS AVERAGE 


CENTERLINE 


———— ARCFLOI 
—O— SURVEY DATA 


—Cl}— #$£=ARCFLOII 
—O— SURVEY DATA 





0.0 0.2 0.4 0.6 0.8 1.0 ~~ 0.0 0.2 0.4 0.6 0.8 1.0 
r/R r/R 


(c) Enthalpy Profile (d) Energy Flux Profile 


Figure 12 Comparison of ARCFLO II and survey data profiles of 
8-cm constricted arc. 
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